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We report about the energy and momentum dependent charge response in Cai.gNao.iCuChC^ 
employing electron energy-loss spectroscopy. Along the diagonal of the Brillouin zone (BZ) we find a 
plasmon peak — indicating the presence of metallic states in this momentum region — which emerges 
as a consequence of substantial spectral-weight transfer from excitations across the charge-transfer 
(CT) gap and is the two-particle manifestation of the small Fermi pocket or arc observed with 
photoemission in this part of the BZ. In contrast, the spectrum along the [100]-direction is almost 
entirely dominated by CT excitations, reminiscent of the insulating parent compound. We argue 
that the observed polarization dependent shape of the spectrum is suggestive of a breaking of the 
underlying tetragonal lattice symmetry, possibly due to fluctuating nematic order in the charge 
channel. In addition we find the plasmon bandwidth to be suppressed compared to optimally doped 
cuprates. 

PACS numbers: 71.30.+h, 73.20.Mf, 74.72. Jt, 79.20.Uv 



I. INTRODUCTION 

Cuprates are still the paradigm for many topics in 
condensed matter physics that cannot be reconciled within 
the conventional Fermi liquid approach. At the heart of 
the research on their electronic structure lies the question 
how the insulating parent compounds develop into high- 
temperature superconductors and finally into ordinary 
metals upon doping. 

It is well established that, upon introducing charges 
into the Cu02 plane cuprates undergo a metal-insulator 
transition (MIT) which is seen in various transport and 
spectroscopic experiments 1 . This MIT, however, appears 
to be rather peculiar because additional charges are be- 
lieved to aggregate in particular patterns termed stripes 
or checkerboards 2 , a phenomenon that is also observed 
in other transition-metal systems like the nickelates 3 . In 
general, the charge order is accompanied by a correspond- 
ing structure in the spin sector . In this respect the very 
recent reports about the existence of charge-only ordering 
in i?Ba2Cu306+a; (R =Y,Nd) 5,6 are remarkable. 

The Ca2- a; Na :r Cu02Cl2 system exhibits charge order 
in the form of bond-centered stripes or nematic behav- 
ior as revealed by scanning-tunneling microscopy (STM) 
experiments 7-9 . Naturally, the role of those inhomo- 
geneities for the superconductivity or more generally for 
the electronic properties of underdoped cuprates is un- 
der strong debate. In particular it is known that in the 
La-based families there is a significant suppression of the 
superconducting transition temperature for x = 0.125 
where the stripe order is most robust 4 ' 10 . While the spin 



dynamics in this region of the phase diagram is well es- 
tablished by inelastic neutron scattering (see e.g. 11 ) the 
corresponding behavior for excitations in the charge chan- 
nel as a function of energy and momentum is, though 
of urgent interest for a complete characterization of the 
collective-mode spectrum, only scarcely explored so far. 



Inelastic electron scattering also known as electron 
energy-loss spectroscopy (EELS) in transmission is a well 
established and bulk-sensitive tool to investigate the en- 
ergy and momentum dependence of the charge response 
in solids as its cross-section is directly proportional to 
Im( — l/e(o>, q)) — with e(uj, q) the complex dielectric func- 
tion of the sample — and therefore allows to investigate 
the dynamics of collective charge excitations 12 . 



Here we employ EELS to study the charge dynamics 
in the system Cai.gNao.iCuC^C^. We find an angular 
dependent transfer of spectral weight between the exci- 
tations across the charge-transfer (CT) gap and those 
associated with the metallic state. The result of this is 
a "nodal metal" that is observed as the truncated Fermi 
surface (FS) arc or small pocket in angle-resolved pho- 
toemission (ARPES) on underdoped cuprates. As the 
observed intensity modulations of the spectrum are not 
compatible with the underlying tetragonal lattice symme- 
try we argue that this spectral-weight transfer appears to 
be accompanied by nematic order in the charge channel. 
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II. EXPERIMENTS AND RESULTS 

The single crystals of Cai.gNao.iCuO^C^ obtained by 
a flux method under high pressure 13 were cut into thin 
id ~ 100 nm) films using an ultramicrotome under ni- 
trogen atmosphere to avoid sample damage due to air- 
exposure. In the spectrometer the films were aligned in 
situ with electron diffraction showing the high quality of 
our samples and allowing for polarization dependent in- 
vestigations along well defined directions within the CuC>2 
plane. The measurements were carried out using a ded- 
icated transmission electron energy-loss spectrometer 14 
equipped with a helium flow cryostat employing a primary 
electron energy of 172 keV and energy and momentum 
resolutions of AE = 80 meV and Aq = 0.035 A~ , re- 
spectively. 

In Fig. 1 we present the behavior of the EELS intensity 
for Cai.gNao.iCu02Cl2 in two high-symmetry directions 
within the CuC>2 plane. The momentum evolution in the 
region of the CT peak around 2 eV were investigated 
in a previous report 15 and we therefore focus on the en- 
ergy range around 1 eV in the following. As can be 
seen for momentum transfers parallel to the [110] direc- 
tion there is a well-pronounced peak in the tail of the 
zero-loss line which disperses to higher energies upon 
leaving the center of the Brillouin zone (BZ). This peak 
is, however, strongly suppressed for g|| [100] (along the 
copper-oxygen bonds). Note that this anisotropy is ab- 
sent in the insulating Sr2Cu02Cl2 16 , in the strongly un- 
derdoped Ca! gsNao.osCuO^CLj 15 and also the optimally 
doped Bi 2 Sr2CaCu208+«5 17 . 

To further quantify the 1 eV feature we present its 
momentum evolution in Fig. 2 for both lattice directions. 
In order to obtain an unbiased estimate for the energetic 
position we did not remove the quasi-elastic line. While 
the energy values for q\\ [110] shown in Fig. 2 track the 
position of the peak maximum around 1 eV this procedure 
could not be applied to the spectra in the [100] direction 
as there the peak is hardly observable. Therefore we took 
the zero crossing of the second derivative between 1 eV 
and 1.5 eV as the characteristic feature for the spectra 
and present its momentum dependence in Fig. 2. Conse- 
quently the different onset energies may be considered as 
an artifact of the data evaluation. As can be seen from 
this analysis, in both directions the dispersion is positive 
and has a bandwidth of about 200 meV. 

To gain a deeper understanding we measured the EELS 
intensity for different temperatures and angles within the 
Cu02 plane. A summary of these results is presented in 
Fig. 3. The left panel shows an angular map measured 
over one quadrant of the BZ at room temperature with a 
constant momentum transfer of q = 0.1 A . Obviously, 
there is a significant angular range where the 1 eV feature 
gains substantial weight at the expense of the CT peak. 
Though not perfectly symmetric 18 , the portion of the BZ 
where the low-energy peak is most pronounced clearly cor- 
responds to the so called nodal direction; the momentum 



region where the superconducting and also the pseudo- 
gap (PC) are known to approach zero from ARPES 19 . 
The temperature (in) dependence of the discussed effect 
is summarized in the right panel of Fig. 3. Those curves 
have been obtained by measuring angular cuts at constant 
energy- loss of 1 eV and 2.4 eV. Subsequently, these two 
spectra were normalized by a similar one measured at 
4 eV energy-loss. The angular independence of the EELS 
signal on the high-energy side allows us to consider this 
third spectrum as a suitable normalization background. 
Note, however, that the observed asymmetric signal does 
not depend on the energy range taken for the normal- 
ization. Clearly, the obtained pattern shows a fourfold 
symmetry over the entire angular range of the BZ with 
the maximum of the 1 eV peak always located around the 
nodal directions (indicated by the vertical dashed bars), 
in agreement with the left panel of Fig. 3. From this 
we conclude that the observed effect is intrinsic and not 
caused by the sample preparation. Were this the case 
the preparation procedure would produce a distinguished 
axis parallel to the cutting direction of the sample, the 
result being a pattern with a periodicity of n instead 
of the 7r/2-period we observe. In addition, the intensity 
redistribution between the features around 1 eV and 2 eV 
is surprisingly robust against temperature variations as 
we do not observe any changes of the periodicity or the 
amplitude between the highest and lowest measured tem- 
peratures shown in Fig. 3 (at least with the resolution 
accessible to us). The same holds true for all intermediate 
temperature steps investigated. 

We emphasize that we do not find evidence for super- 
structure reflections in the elastic scattering neither at 
high nor at low temperatures. This agrees with earlier 
x-ray scattering experiments on this system 20 . In addi- 
tion the pattern seen in Fig. 3 changes quantitatively (see 
below) but not qualitatively when changing the momen- 
tum transfer between q = 0.08 A and q ~ 0.4 A 
and we do not see any indication for a resonance-like 
enhancement of the effect at the momentum transfer of 
the checkerboard q = 2n/AaQ ps 0.41 A 7 (see Fig. 4) 
with the lattice constant do = 3.85 A. 

In order to track the momentum evolution of the inten- 
sity anisotropy we plot the spectra along the [100]- and 
[110]-direction for a series of q- values in Fig. 4. Already 
by visual inspection it is clear that the anisotropy persists 
but is not monotonic function of q. 

As a more quantitative measure of the anisotropic EELS 
intensity we define the following ratio: 

r>SW _ bW <l ( 1 \ 

q swl wo] 

where the "spectral weights" SW* are determined by 
SW* = f"* duiu) I{oj,q) with the normalized intensities 
I(u>, q) shown in Fig. 4 and the momentum transfer q 
parallel to the directions indicated by the superscript 
i € {[100] , [110]}. Note that if one converted the in- 
tensities I(u>,q) to the absolute value of Im(— l/e(w, q)) 
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Figure 1. (Color Online) The EELS intensity in Cai.gNao.iCuC^Cb for momentum transfers parallel to [100] (left panel) and 
[110] (right panel) measured at room temperature. Note the intensity enhancement around 1 eV parallel to the [110] direction. 
All spectra have been normalized on the high-energy side between 3.5 eV and 4 eV. 
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Figure 2. (Color Online) The momentum dependence of the 
1 eV feature seen in Fig. 1. See text for details. 



by means of a Kramers-Kronig (KK) transform the in- 
tegrals SWg would indeed correspond to the spectral 
weight (and also be constrained by the /-sum rule 12 ) 
which motivates our nomenclature. We did not, however, 
perform the KK-calculation, because, in particular for the 
[100]-direction, the subtraction of the quasi-elastic line 
is highly ambiguous. The boundaries of the integration 
lo\ and ix>i correspond to the dips in the EELS intensity 
below and above the 1 eV feature along [110], respec- 
tively (see the vertical arrows in Fig. 4). The evolution 
of the ratio defined in Eq. 1 as a function of momentum 
transfer is visualized in Fig. 5. As can be seen, there is a 
rather well defined maximum around the incommensurate 
g c ~ 0.175 A" 1 ... 0.2 A~\ 



III. DISCUSSION 

Although there is experimental evidence that the in- 
frared optical response in underdoped cuprates contains 



more than a simple Drude term 21,22 we attribute the in- 
tensity increase around 1 eV and in particular the well 
pronounced peak in the [110] direction to the Drude plas- 
mon caused by the doped charge carriers. This is moti- 
vated by the metallic resistivity and the clearly visible 
plasma edge seen in the reflectivity of Cai.gNao.iCu02Cl2 
at this energy 21 . In addition earlier EELS reports on op- 
timally doped cuprates, in particular from the Bi- and Y- 
families, found similar features that have been interpreted 
analogously 17 ' 23 ~ 26 . In contrast to our observation from 
Fig. 2, however, all these experiments show a quadratic 
plasmon dispersion that follows the conventional behavior 
predicted for simple metals 27 with a bandwidth that is 
larger than the one reported here by at least a factor of 
two. As the slope of the plasmon dispersion in a sim- 
ple metal is proportional to the squared Fermi velocity 
vp ~ n 1 / 3 with the charge-carrier density n the smaller 
dispersion in Cai.gNao.iCu02Cl2 might simply reflect the 
(compared to the optimally-doped systems) smaller FS 
seen in ARPES 28 . This is also in line with the systematic 
blue-shift of the plasma edge as a function of doping in 
Ca 2 - :c Na x Cu0 2 Cl 2 21 . 

We now turn our attention to the obvious anisotropy 
in the charge response for different directions within the 
Cu0 2 plane. 

According to earlier reports 13 ' 21 ' 29 ' 30 but also from our 
elastic scattering data (not shown) the crystal structure 
of the C^-rrNa^CuC^C^ system is perfectly tetragonal 
without any signs of buckling or orthorhombicity often 
found in other cuprate families. This implies that for 
q = the dielectric function — which is in the general 
case of anisotropic systems a tensor — is diagonal and 
contains only one independent in-plane component and 
consequently should be isotropic throughout the entire 
Cu02 plane 31 . Therefore also l/e(w, q), which we probe 
in EELS, should be describable by a single isotropic com- 
ponent. Although we cannot measure at q = due to the 
increasing influence of the zero-loss peak and the enhance- 
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Figure 3. (Color Online) Left panel: Angular (polarization) dependence of the EELS intensity for q = 0.1 A 1 measured at room 
temperature normalized on the high-energy side between 4 eV and 5 eV. Right panel: Constant energy cuts for q = 0.1 A 1 
around the full BZ (see text for details). Note the fourfold symmetry. The dashed bars correspond to the nodal directions as 
derived from the crystal structure measured with elastic scattering. In both panels the angle is measured relative to the [100] 
direction. 




Figure 4. (Color Online) The asymmetry of the EELS intensity for the two in-plane directions as a function of momentum. For 
each momentum value, the curves are normalized between 3.5 eV and 4 eV and the arrows indicate the integration ranges for 
the evaluation of Eq. 1. 



merit of surface scattering 27 , in general the momentum 
transfer of q — 0.1 A , which corresponds to only about 
6% of the BZ size, is small enough to be considered as 
the optical limit (q — 0) (see also the isotropy for this 
values of momentum transfer in the response of other 
cuprates 15-17 ). Consequently, we attribute the behavior 



described above to an intrinsic (electronic) breaking of 
the underlying tetragonal lattice symmetry as our data 
imply strong polarization-dependent changes of the charge 
response, equivalent to two in-plane components of the 
dielectric tensor even in the limit q — > 0. This is com- 
patible with the STM data 8 ' 9 where evidence is found 
for a local breaking of the symmetry from C4 (tetrago- 
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Figure 5. The momentum dependence of the ratio defined in 
Eq. 1. 



nal) to C2 (orthorhombic) which can be ascribed to the 
presence of a nematic order parameter 2,32 . Note, how- 
ever, that while the STM data implies a higher metallic 
character parallel to the bonds our data seem to indi- 
cate an enhanced density of states along the diagonals 
of the unit cell (BZ). Our data therefore also appears 
at odds with results based on Raman- 33,34 and neutron- 
scattering 35 where stripes are along the diagonal only for 
a doping too low to induce superconductivity. There are 
several possible scenarios to reconcile the observed ir/2- 
periodic signal (see Fig. 3). Either there are perpendicular 
ID domains within a single CUO2 plane or the charge 
modulation is 2D. Another possibility is that different 
CUO2 planes along the path of the electron beam have 
perpendicular ID domains and the observed intensity 
pattern is just a superposition of the contribution of all 
unit-cells in the sample. In any case, along the domains 
there is metallic transport reflected by the appearance 
of the plasmon whereas perpendicular to them the spec- 
trum is dominated by the CT excitations reminiscent of 
the insulator. These domains are highly fluctuating but 
the timescale of the electron-scattering process is short 
enough that EELS measures a snapshot of this dynamic 
behavior. At low enough temperatures these charge fluc- 
tuations may lock in to form a well-ordered pattern but 
due to the transmission geometry of the experiment we 
are limited to T ^ 20 K. From the absence of super- 
structure reflections and the robustness of the tetragonal 
lattice symmetry, we conjecture that Cai.gNao.iCuC^C^ 
tends to show nematic order. Our data therefore also 
indicate that the regular checkerboard pattern seen in 
the STM 7 may be pinned by the surface, in agreement 
with the findings of Ref. 36. From Fig. 5 we identify a 
characteristic momentum q c ~ 0.175 A ... 0.2 A of 
the nematic fluctuations which translates to a length scale 
of l c = — -30A...35A-8a ...9a . 

From the ARPES literature on underdoped 
Ca2_a;Na ;r Cu02Ci2 28,37 and also other cuprate systems 19 



it is known that the large FS observed in the optimally 
and overdoped compounds is confined to a small 
angular range — either in the form of an arc or closed 
pocket 38 — around the nodal direction which increases 
with doping. From our data, in particular from Fig. 3 
it is clear where these metallic states derive from: they 
are created by a transfer of spectral weight from the 
upper Hubbard band (the excitations around 2.5 eV in 
Fig. 3) down to low energies in the vicinity of the Fermi 
level. In EELS these low-lying states show up at finite 
energies (around 1 eV) because for a metal a peak in 
Im(— 1/e) = £2/(61 + e|) occurs at the plasma frequency 
which is determined by the density of free charges and 
therefore finite. Note that such a spectral- weight transfer 
is considered to be a hallmark of strongly-correlated 
systems . Angle-resolved inverse photoemission could 
reveal the corresponding behavior in the single-particle 
channel. 

Interestingly, a similar anisotropic appearance of metal- 
lic states upon doping is also predicted theoretically for 
the Hubbard model by cluster extensions of dynamical 
mean-field theory 40 and variational cluster approaches . 

As the spectral- weight transfer, however, is expected 
to occur symmetrically for all four nodal points of the 
BZ and therefore preserves the underlying tetragonal 
symmetry of the lattice our observations seem to require 
the collaborative action of (fluctuating) charge order and 
spectral- weight transfer. 

IV. SUMMARY 

To summarize we investigated the energy and momen- 
tum dependent charge response in Cai.gNao.iCu02Cl2 by 
means of electron energy-loss spectroscopy. We find a po- 
larization dependent transfer of spectral weight between 
the CT excitations — dominating the spectrum around 
(tt, 0) — and the metallic plasmon (around (tt, it)) which 
can be understood as an anisotropic metal-insulator tran- 
sition within the BZ. The observed pattern of the charge 
response indicates the breaking of the underlying tetrag- 
onal lattice symmetry which we take as indicative for a 
nematic order parameter in the Cu02 plane. In addition, 
the plasmon is found to disperse with a bandwidth smaller 
than in optimally doped cuprates. All effects are robust 
against temperature variation arguing in favor of a highly 
fluctuating order. 
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